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ON , Abstract. The existence of a large number of asymmetric, rotating neutron stars, 

each individually emitting periodic or quasi-periodic gravitational waves in the fre- 
£h , quency band around 100 Hz, raises the possibility of detecting their combined signals, 

by exploiting the amplitude modulation of the received waves as the antenna changes 
its orientation with respect to fixed stars. This modulation is directly related to the 
amount of anisotropy present in the source distribution, and, if detected, could give 
valuable information about the spatial distribution of neutron stars in our Galaxy. 
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Introduction 



It is possible that a fraction of the gravitational wave (GW) stochastic back- 
' ground is generated astrophysically within our Galaxy, e.g., by a large popula- 

tion of galactic sources emitting long-duration individual signals which overlap 
in frequency. This is the case, for example, in the low-frequency region, where 
the expected signal from unresolved close binary systems represents an effective 
Oh' sensitivity limit for planned spaceborn interferometers, at frequencies below ~ 3 

mHz|J. 

S-h_ | A similar situation may arise for ground-based interferometers, when con- 

. sidering waves emitted by a large population of non-axisymmetric, rotating 

neutron stars (NS). This contribution to the background signal, in any case, 
would not have the same relevance as the galactic binaries background has for 
' space interferometers, for the following reasons: first, we do not know how many 

• rotating NS are present in our Galaxy, nor how efficiently they emit GW, and 

' therefore we do not know whether their contribution to the background sig- 

nal is relevant, when compared to other stochastic sources. Also, one should 
consider that the stochastic nature of this signal, in the sense that individual 
sources cannot be resolved, results from the finite duration of the observation 
time (i.e. from the finite size of the resolution bin in the frequency spectrum), 
which in the terrestrial case does not have any a priori upper limit. Finally, 
even if one assumes that the number of sources is large enough for the emitted 
waves to be unresolved in a given integration time, one can always rely on a 
network of detectors, which would allow disentangling this signal from other 
stochastic sources. One can thus conclude that the overall signal from galactic 

t Paper presented at the 2nd Workshop on Gravitational Wave Data Analysis, Orsay 
(France), 13-15 November 1997. 
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NS should not interfere with the search of other sources, unless the NS birthrate 
is considerably higher than today's estimates (see Eq. (||) below). 

It has been suggested, however, that detecting and measuring this signal 
with a single detector will provide direct and valuable information about the 
spatial distribution of NS in the Galaxy_B, and in some cases could even be 
easier than detecting individual sources || . This detection strategy is based on 
the fact that the antenna response is not isotropic, and therefore the change 
of orientation of the interferometer with respect to fixed stars generates an 
amplitude modulation of the stochastic signal, given that the sources are not 
distributed isotropically with respect to us. The important point is that this 
modulation is totally deterministic, and for a terrestrial interferometer it can 
be easily expressed as the sum of four sinusoidal oscillations, with frequencies 
k/Td (k = 1, . . . ,4), where Td is the duration of the sidereal day The ampli- 
tude of these four periodic oscillations is related to the amount of anisotropy in 
the spatial distribution of sources, being identically zero for a perfectly isotropic 
population (as it should be, since in this case the detector's motion is totally 
irrelevant). Thus, by measuring the modulations of the variance of the stochas- 
tic signal one should be able to obtain valuable information about how the NS 
are distributed in our Galaxy g. 

This contribution summarizes some recent results on this subject. Full de- 
tails of the calculations can be found in la] . 



2 The strength of the stochastic signal 

In general, deformations produced by internal forces in a slowly-rotating NS are 
responsible for creating a slightly non-axisymmetric configuration (described by 
the equatorial oblateness e e ), and for making the star's angular velocity precess 
around the near-symmetry axis (we call 'wobble angle' 9w the constant angle 
between the rotation and symmetry axes). Both effects result in the emission of 
sinusoidal gravitational radiation, respectively at twice the rotation frequency 
/i = uj rot /TT, and at its precessional sideband f% = (cj rot + uj prec ) / (2n) ~ /i/2. 
The amplitudes of the two components are || 

32 2 e e Ifi 

h m t = y -jir , (1) 

[32 2 e p 8 w I{2f 2 f 

Vcc = y y7r , (2) 

where / is the star's moment of inertia, e p is the poloidal oblateness, mainly 
due to centrifugal forces, and r is the distance to the earth. 

Each NS is slowly spinning down due to the loss of energy and angular 
momentum in the form of electromagnetic and gravitational waves. The exact 
evolution of the star's rotation frequency, and consequently of the GW frequency 
and amplitude, depends on which of the two emission mechanisms is dominant. 
This also affects in a significant way the spectral properties of the cumulative 
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signal from the whole population of NS. For example, for a given birthrate v and 
a given observing time T b s , there exists a threshold frequency f c below which 
the signal from the NS population behaves like a stochastic signal, in the sense 
that one expects, on average, more than one source per frequency resolution bin. 
It turns out that the relatively fast evolution due to electromagnetic emission 
makes this continuum frequency f c being located below the frequency band of 
interest (/ c < /,„;„ ~ 20 Hz), unless 

v >2x 10 3 yr"\ (3) 

a birthrate considerably higher than the present estimate v es t — 10~ 2 yr~ 1 . 

If, however, a significant fraction of sources are born with a small magnetic 
field, so that GW radiation reaction dominates the (much slower) evolution, 
then the limit (^) drops to 

v >5.5 x lO^V^ 1 ■ (4) 

In deriving eqs. @ and (|J) it has been assumed a magnetic field strength 
\B\ ~ 10 12 Gauss, an ellipticity e ~ 10~ 6 , a moment of inertia / ~ 10 45 g cm 2 , 
and T b s = 1 yr; a more general expression for the minimum birthrate v can be 
found in ||. 

In view of the fact that unresolved NS are likely to be spinning down by 
gravitational radiation reaction, we will consider only this mechanism from now 
on. Also, for the sake of brevity, the precessional component /i prC c of the radi- 
ation will be neglected, by imposing 9w = 0. The rotational component h TO t 
gives, in the confusion region / < f c , a stochastic signal of spectral amplitude 



S h (f) = 77 / ^rdr ~ 1.5 x 10" 51 Hz 



5/ J r 2 V 100Hz 



v 



10-2 yr" 1 / V10 45 g cm 2 



J r z 

V 

(lOkpc)" 2 



(5) 



Note that the spectral level (^) does not depend on the oblateness of the 
sources, a direct consequence of the fact that gravitational radiation reaction is 
assumed to be the dominant spindown mechanism. 

As one can infer from Eq. (||), a birthrate of one pulsar each 100 years gives 
a signal well below the expected sensitivity of planned terrestrial detectors. Un- 
fortunately, the actual rate is rather unknown, since a few hundred pulsars out 
of perhaps 10 8 rotating NS in our Galaxy have been detected, and presumably 
the strongest sources have not yet been detected electromagnetically. However, 
as shown in the next section, this stochastic signal can be detected even if its 
strength is much less than that of the instrumental noise, by exploiting the 
peculiar time dependence of the signal's amplitude. 
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3 Amplitude modulation of the stochastic signal 



As already mentioned in the Introduction, the stochastic signal generated by 
any anisotropic distribution of sources, as seen by an interferometer anchored 
on the surface of the earth, is non-stationary. It is now possible to quantify this 
effect, and relate it to the actual distribution of sources via a simple integral. 

In particular, we will focus our attention to the signal's expected variance, 
which can be expressed as 

a 2 (t) = JdfS(f,t), (6) 

where S(f, t) is the generalization of Eq. (^|) which takes into account the actual 
response of the interferometer, through its antenna pattern P(0, cj>, i), i.e. 



2/ 

v 

Note that Eq. (|t]) becomes identical to Eq. (jjj) if the spatial distribution is 
isotropic with respect to us, since JdCl P(8,(/),t) = 8tt/5. For a terrestrial 
interferometer in the equatorial reference frame one finds 

4 

P(0, 4>,t)=J2 { a k{0, 0) cos{kLu d t) + b k {9, (f>) sm(kcu d t)} , (8) 

fc=0 

where the Fourier coefficients a k , b k are explicitly given in B , and Lo d is the 
earth sidereal frequency u> d — 7 x 10~ 5 rad s _1 . Inserting Eq. in Eq. (Q) one 
trivially obtains, upon integration over frequency 

4 

o 2 {t) = ^ i a k cos(kuj d t) + f3 k s'm(kuj d t)} , (9) 

fe=0 

where 

oo 

a k tx Jd£la k (9, <j)) J dr p(r,6,c/>) , 

4-7T 

and analogous for f3 k , with a k replaced by b k . Thus, the variance of the stochas- 
tic signal is a deterministic function of time; its Fourier decomposition reveals 
four discrete frequencies, harmonics of the earth sidereal frequency uj d . The 
amplitudes of these four spectral lines depend on the particular form of p, i.e. 
on the spatial distribution of the sources. By measuring them, one can therefore 
obtain information about the global distribution of sources in the Galaxy. 

The detection strategy could be as follows: the data record, of duration 
T bs 3> T d , will be divided into several batches of duration r <C Tj. For each 
batch J we obtain an estimate Sj of the variance a 2 (tj). According to Eq. (^), 
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the variables {sj} are expected to show four sinusoidal modulations, of known 
frequencies (ujd and its first three harmonics) and phases (depending on the 
initial orientation of the detector with respect to the Galaxy). The amplitude 
of these oscillations can be estimated with a simple Fourier analysis, in the 
hypothesis that the instrumental noise is stationary. Note that, in order for 
this search to be successful, it is not required that Sh > S no i se , since we are 
not directly comparing the GW stochastic signal against the instrumental noise; 
instead, we are trying to detect a deterministic and periodic signal a 2 (t) buried 
in a stationary random noise of amplitude cr 2 loise . The main problem in this 
measurement comes from those noise sources that, far from being stationary, 
present the same diurnal variations as the gravitational signal, due for example 
to temperature changes in the mirrors. This issue has been addressed in great 
detail by Giazotto et al. where a careful monitoring of the temperature 
fluctuations in the mirrors has been advocated. 

In order to estimate the magnitude of the effect, and in particular the am- 
plitude of the periodic terms of a 2 (t) with respect to the constant one, we have 
considered different galactic populations, like the bulge, the disk, the halo, etc. 
Table 1 shows some results in this direction: for each of the three interferome- 
ters presently under construction - the two LIGO detectors in the United States 
[Q, and the VIRGO detector in Italy -, we have estimated the amplitudes 
7^ = y/ a\ + Pi (k — 0, . . . , 4) generated by, respectively, (i) a spherical halo, 
(ii) an exponential disk, and (iii) a cluster of sources located in the galactic 
center. 
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Table 1: Summary of the Fourier coefficients 7fc = \J o? k + /3| , for various galac- 
tic components. For each population, the columns marked L-l, L-2, and V are 
obtained taking into account the orientation on the earth surface of, respec- 
tively, the two LIGO detectors and the VIRGO detector. For further details, 
see [||. 

From Table 1, one can conclude that the largest contribution to the periodic 
terms of a 1 comes from the galactic disk. Also important is the signal coming 
from a localized cluster of sources, like the galactic bulge or even a rich globular 
cluster. The periodic contributions from the spherical halo are expected to be 
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negligible. 
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